Implications of the Intermediate Mass Black Hole in globular cluster Gl on Dark 

Matter detection 
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Recently there has been a growing evidence in favor of the presence of an Intermediate Mass Black 
Hole in the globular cluster Gl, in Andromeda Galaxy. Under the assumption that formation of 
this globular cluster occurred within a dark matter halo, we explore whether the presence of a black 
hole could result in an observable gamma ray signal due to dark matter annihilation in this globular 
cluster. Starting from an initial NFW matter profile, with density parameters consistent with Gl 
observations, we find that indeed, if the spike in the density has been formed and has survived till 
present, the signal could be observed by GLAST and current ACT detectors. 



PACS numbers: 95.35.+d 



I. INTRODUCTION 



One of the promising ways of detecting self- 
annihilating dark matter is through its final annihila- 
tion products, gamma rays being one of them. These, 
so called indirect detection techniques, offer an oppor- 
tunity to identify dark matter (DM) properties, comple- 
mentary to other searches, as dark matter direct detec- 
tion and accelerator experiments. Gamma rays with en- 
ergies higher than 1 GeV have been searched for, in both 
satellite based (EGRET [lh and ground-based experi- 
ments (HESS [|, MAGIC Hi, CANGAROO Q, VERI- 
TAS [5j,...) and new telescopes are being planned and 
built, such as GLAST [| with launch scheduled for the 
next spring. The gamma ray signal of self-annihilating 
dark matter depends on the density squared and there- 
fore "hot spots" where DM density is enhanced, are ex- 
pected to exist across the sky. Promising sites for gamma 
ray DM detection should, therefore, have high DM den- 
sity, be located in our vicinity and have a low level of 
astrophysical background, as the centers of our own and 
nearby galaxies, dwarf galaxies of the Milky Way, earth 
sized halos or primordial intermediate-mass black holes, 
for a review see Q ■ 

In this paper we propose yet another site for dark mat- 
ter indirect observation, the globular cluster Gl in the 
Andromeda Galaxy. This cluster has recently been ob- 
served with X-ray Q and radio @ telescopes and both 
measurements detected a source in the cluster core, con- 
sistent with the presence of an Intermediate Mass Black 
Hole (IMBH). Dynamical considerations, based on the 
most recent HST/WFPC2 and KECK/HIRES spectro- 
scopic measurements support an IMBH presence, [To[ 
(see also, [n], [lU, HH). In this letter we explore whether 
the presence of an IMBH could have an effect on the DM 
density favorable enough for Gl to be a promising target 
for indirect detection. Dark matter halos react to the 
growth of black holes. In the case of adiabatic growth 
of a central black hole, the density is strongly enhanced 
in its vicinity, resulting in a spike [Til ]. Since the gamma 
ray flux depends on the square of the density, a spike 
would enhance the flux by several orders of magnitude. 
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The existence and the slope of a spike is highly uncertain, 
however. We will later see how spikes with quite shallow 
slopes could still produce an observable signal. 

Gl as a dark matter observation target suffers from 
being rather far from the Earth, but as we will see below, 
the enhancement of the density due to a black hole can 
potentially compensate for the distance. For example, 
if compared to the dwarf galaxy Draco, Gl is 10 times 
farther, but the signal from Gl, in the case of a spike 
could be 10 2 ^ 4 times stronger than the one expected from 
Draco, 15j. Also, on the positive side, globular clusters 
are not expected to have significant background in high- 
energy gamma rays, as we will argue later. Relic IMBH 
have long been considered good sites for DM detection, 
but suffered from uncertainty in their location. If Gl 
really hosts an IMBH it would offer a unique opportunity 
to study these objects and their interaction with dark 
matter. 



II. GLOBULAR CLUSTERS AND Gl 

Globular Clusters are very dense stellar systems: they 
have the typical mass of dwarf galaxies, but the size a fac- 
tor of ~ 10 smaller. They are found in the halo or bulge 
regions of galaxies, Because they are old, they do 

not have significant amounts of interstellar gas and the 
stars within globular cluster are usually coeval. This sug- 
gests that the astrophysical background in high energy 
gamma ray telescopes should be small. 

In the primordial scenario of globular cluster formation 
[Tl| . globular clusters formed in DM minihaloes, before 
or during reionization, being the first galaxies to form; It 
is plausible, therefore, that globulars initially had a dom- 
inant dark matter component. They were then accreted 
by larger galaxies and lost their extra nuclear material 
due to stripping by the tidal forces. Today globular clus- 
ters are baryon dominated, with M/L of a few. They 
have very dense stellar cores and the presence of tidal 
tails observed on globular cluster Palomar 5 suggests that 
they do not have an extended DM halo [l8|. How much 
dark matter could be present in globular clusters? King- 
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Michie models do a satisfactory job in the dynamical 
description of globular clusters while considering only its 
stellar population. In that sense, presence of exotic dark 
matter is not required and no lower bound on DM an- 
nihilation signal from Gl can be placed. However, the 
simulations of Mashchenko and Sills [2(| of the evolution 
of globular clusters in tidal fields show that the dark mat- 
ter present in their nucleus could be compressed and its 
density profile steepened by these dynamical processes 
and that, depending on the initial dark matter profile, 
and type of DM-stellar interaction, some globular clus- 
ters could still present an extended halo. 

The Gl cluster is one of the most luminous and the 
most massive globulars in the Local Group. Its total lu- 
minosity is L ~ 2 10 6 Lq and its mass is estimated to 
be between 1.4 and 1.7 x 10 7 M Q in King-Michie models, 
PH . However, the values of the total mass of the cluster 
has large scatter depending on the model used for its cal- 
culation. For example Baumgardt et al. (l3j find a value 
of (8±1) 10 6 Mq in their simulation which uses an evolu- 
tionary model. The main reason for these discrepancies is 
that spectroscopic data on Gl are rather poor: measure- 
ments of the velocity dispersion profile in the cluster was 
performed only in the inner region of the cluster (within 
*~ 5 pc, [l(|) or the velocity dispersion was reduced to a 
single value, measured within a slit ~ 6 x 24 pc, [22J. 

Gl has a stellar core radius of 0.52 pc, and a half mass 
radius of 13 pc; its tidal radius is about 200 pc, [2l|. The 
cluster lies at a distance of 770 kpc from us and it would 
be clearly separated from the bright Andromeda's disk 
in current high energy gamma ray telescopes. 

Gl has several properties not typical for globular clus- 
ters. It has a very large central velocity dispersion [22T ]. 
has a very flattened shape, with mean ellipticity of 0.2 
and it shows a spread in metal abundances among its 
stellar population. In fact, one of the possible scenarios 
of Gl formation is that it is a surviving nucleus of a dwarf 
elliptical galaxy, which would have lost its outer envelop 
through tidal interaction with M31. In [23|, Bekki and 
Chiba show that also in this scenario presence of DM in 
the core is plausible. 

Evidence for a black hole in Gl: While there is a grow- 
ing knowledge of super-massive and stellar mass black 
holes, intermediate mass black holes are yet poorly un- 
derstood. Recent detections of AGNs in low-luminosity, 
late-type galaxies suggest that IMBHs do exist (see, for 
instance, [24|). 

Several papers dealt with a dynamical analysis of the 
Gl cluster and the possib le p resence of an IMBH with 
different outcomes fiol. Jl2l [ilij . The most recent analysis 
by Gebhardt et al. [1J| (G05), favors its presence. They 
find the best fit black hole mass to be 1.8 ±0.5 x 1O 4 M . 

X-ray observations with XMM Newton reported a 
source with luminosity expected from an accretion onto 
an IMBH [25|. Radio observations, using Very Large 
Array detected a faint radio source within an arcsecond 
of the cluster core [j| with luminosity consistent with 
mass estimates by G05 and measured X-ray luminosity 
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FIG. 1: The non-shaded region is the range of r s and mnu 
which could be detectable by GLAST or VERITAS with 5a 
significance or better. Shaded region is the excluded region 
from the EGRET observation. Solid line represents the case 
of AM = 9 10 4 M Q , and dotted one AM = 3 10 4 M©. 



thereby strengthening the evidence for a black hole pres- 



III. DARK MATTER DENSITY PROFILE 

We consider a cuspy NFW dark matter profile, p(r) = 
p$(r / r 8 )~ l (\ + r/r s )~ 2 where po is central density and 
r s the scale radius. 

Dark matter halos react to the growth of black holes, 
and in the case of adiabatic growth, the result is the 
formation of large over-densities called spikes, If the 
initial dark matter density was cuspy, after the adiabatic 
growth, the dark matter assumes a distribution given by 



Ps P ,NFw{r) = Pcus P (R sp ) g~t{r) ( — ^ 
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(1) 



where R sp is the radius of influence of the black hole, 
found as a solution of the equation J^ 3V p(r)d 3 r = 2Mbh 
g 1 (r) ~ (1 — 4Rs/r) 3 , where Rs is the Schwarzchild ra- 
dius of the black hole and 7 sp = 7/3 in the NFW case. 

On the other hand, in models with finite cores, the 
slope of the spike is shallower, j sp = 3/2. Numeri- 
cal simulations of halo formation predict cuspy profiles, 
0. However, simulations do not take into account the 
baryon-DM interactions and have other limitations, and 
we therefore consider a Burkert profile [26|, which has a 
core at small r, p(r) = po(l + r/r s )~ 1 {l + (r/r s ) 2 ] -1 . We 
parametrize the final dark matter profile after the growth 
of the spike as p sp .core{r) = p , cor e 9j(r)(R sp /r) 3/2 . 

Due to annihilations in the inner regions of the spike, 
there will exist a maximal dark matter density, given 
by Pmax = Tn, I {o~v)tBH , where m is dark matter mass, 
(o~v) is the thermally averaged annihilation cross section 
and tBH is the age of IMBH, which we take to be 10 10 
years. The final dark matter profile is then given by 

— Psp Pmax i \Psp H~ Pmax)- 
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The adiabatic growth approach used here maximizes 
the effects of the black hole on the dark matter density 
profile. Substantial work has been done on a more realis- 
tic description of black hole influence on dark matter evo- 
lution, and the existence of a spike in our Galactic Center 
has been questioned, [13] • Various processes, such as off- 
center black hole formation, black hole merging events, 
dark matter scattering by stars, etc. could all result in 
slopes shallower than predicted by the adiabatic growth 
model. The reconstruction of the central mass profile of 
Gl, in G05, points to the existence of a cusp, both in the 
central luminosity and in the mass profiles. This in turn 
suggests that Gl did not go through a recent merger, 
which would deplete both stellar and dark matter cusps. 
On the other hand, stellar-dark matter scatterings might 
be important in Gl. Gravitational scattering of dark 
matter with stars could deplete the spike and the slopes 
are generically well described with slope 3/2, if this ef- 
fects are taken into account, [28|. To probe this effect we 
somewhat artificially vary 7 sp and find the lowest value 
which would result in a detectable signal. 

To construct the DM density in Gl we start from 
the assumption that the outer part of the dark matter 
halo has been stripped together with the stellar envelope, 
while the dark matter in the center has survived due to 
the compactness of the core, as motivated by above men- 
tioned simulations, [13, HH. We conservatively truncate 
the profile at the tidal radius. Since the main contri- 
bution to the annihilation signal comes from the inner 
region of the profile {r<r s ) this assumption would not 
significantly affect the result. Also, we conservatively as- 
sume that the inner slope of the density profile does not 
change with the evolution of the cluster. 

The two parameters of the density distribution, r s and 
Po, are poorly constrained in globular clusters. One typ- 
ically constrains r s and po by solving spherical Jeans 
equation, based on the velocity dispersion profile mea- 
surement, [29L |30|. In the case of Gl, G05 find the mass 
density profile in the inner 5pc, using this approach (the 
cllipticity of Gl is less than 0.1 in the innermost region of 
the cluster, [2l[, suggesting that the spherical approxima- 
tion might be valid there), but determination of r s and 
po is further complicated by the fact that Gl is baryon 
dominated, and that the stellar M/L dependence on the 
radius is not well understood. We choose r s and p by 
making sure that the dark matter density profile satisfies 
the three constraints listed below: 

i) We use the uncertainty on the total cluster mass 
AM as the limit on the amount of dark matter allowed 
in the cluster: J pDMcPr = AM. As AM we use both 
the error on the total mass within the King-Michie model 
(3 1O 6 M ) or the mass dispersion between various models 
(9 10 6 M Q ). 

ii) When adiabatic compression due to the black hole 
is present, the strongest limit might come from the mea- 
surement uncertainty of the total mass measured in the 
first bin (~ 0.1 pc), which is - 3 10 3 M Q , [nj. Wc use 
this mass as an upper limit on the amount of DM in the 



central ~ 0.1 pc; 

iii) We use the mass density profile, as found by G05, 
by imposing that the dark matter density should not be 
bigger than 10% of the found mass density of the cluster 
at any distance (10% is somewhat arbitrary value, error 
bars on the mass density of the cluster as calculated in 
G05, are hard to estimate; the measurement error in the 
central bin is 15%). 

For comparison, we notice, that the values of r s (and 
po) can be estimated also by using the virial mass 
and the concentration parameter c V i r , (for details, see 
[31I]). If we take for M OT > to be ~ 5 times the mass of 
Gl today, it leads to a value of r s for the NFW profile of 
~ 270 pc, which is close to its tidal radius today. 

IV. DETECTION SIGNAL 

In this section we estimate the strength of the dark 
matter annihilation signal as expected in the GLAST and 
ground-based Atmospheric Cerenkov Telescopes (ACTs). 
The EGRET experiment did not measure any signal 
above the diffuse gamma ray background in the direc- 
tion of Gl. We use this fact to place further limits on 
the value of r s and po- 
rn the case of Gl, when the typical scale of the dark 
matter halo is much smaller than the distance to the ob- 
ject, the flux of gamma r ays p roduced by self-annihilating 
dark matter is given by, [la ]. 

4nD 2 J 2 \m dm ) K ' 

where D is distance to the object and mdm is the dark 
matter mass; iV 7 is number of photons produced per an- 
nihilation and it is calculated here by means of the ana- 
lytic fitting formula for the 7 ray spectrum produced by 
particles which annihilate dominantly to gauge bosons 
[32j which is a sufficiently accurate approximation for 
our purposes; we assumed a value for DM annihilation 
cross section (av) of 3 10 -26 cm 3 s _1 , corresponding to 
the value expected for thermally decoupled relic. 

For GLAST experiment we assume exposure of 1 m 2 yr. 
The biggest contribution to the background for Gl comes 
from the Galactic diffuse emission at intermediate lat- 
itudes. Its strength is <6 10~ 6 (GcV/E)- 2 cm" 2 s _1 
sr" 1 GeV" 1 , |3]. Using field of view of 10" 5 sr, this 
background accumulates to about 20 events in the given 
exposure. 

For ACT experiments we assume an effective area of 
4 10 s cm 2 , 40 hours of observation and an energy thresh- 
old of 100 GeV. The dominant background in ACTs 
comes from hadronic showers which are misidentified as 
gamma rays. Its spectrum is given as dNhad/dE ~ 
3 ew(£VGcV" 2 - 7 ) GeV- 1 sr" 1 cm" 2 s" 1 , |ljf, where 
£had is the fraction of misidentified hadronic showers. For 
Chad = 0.01, typical for ACTs, and a field of view of 10 -5 
sr, this background accumulates at a rate of ~ 80 events 
per hour of exposure. 
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In the direction of Gl, EGRET observed flux of ~ 
5 10~ 7 cm 2 s _1 sr _1 above 1 GcV, The area of 

EGRET was 6400 cm 2 , the angular resolution 0.5° and 
we will assume one year of effective time. The measured 
number of events, from the direction of Gl, in the energy 
range 1-10 GeV is ~ 25. 

We show the results in a Figure [TJ where a NFW pro- 
file is assumed. The lower shaded region is excluded by 
EGRET. The upper right corner of the figure is observ- 
able by either GLAST or VERITAS experiments, with 
5cr significance, or better. We did not consider values of 
r s higher than 300 pc, since the tidal radius of Gl is 200 
pc. We comment here that spikes with slopes as shal- 
low as 7 S p = 1.75 could still be observable, albeit for a 
very small parameter space: say for m — 100 GeV, the 
observable range of r s would be 3 pc <r s <5.2 pc (for 
AM = 3 10 6 A'/q), but with the parameter space enlarg- 
ing considerably for higher j sp . This is meant to illus- 
trate that even slopes shallower than the ones predicted 
by adiabatic growth of a black hole, could potentially be 
observable. We did not perform a detailed analysis in 
this case, instead we varied ad hoc the values of "f sp in 
Eq. (fT]) in order to mimic complicated superposition of 
effects in the central DM region (initial steepening of DM 
cusp, motivated by [2(| , and various dynamical processes 
that might contribute to the depletion of the spike). If 
the dark matter profile contains no spike, we find that its 
selfannihilation signal would be unobservable by GLAST 
and current ACTs. We also find that the core profile 
would not produce an observable signal even if the adi- 
abatic spike is present. The main reason for this is the 
combined effect of a shallow slope of a spike and strong 
constraints on DM density profile parameters in the case 
of core profile. 



V. SUMMARY 



In this paper we show that the dark matter annihila- 
tion signal from the globular cluster Gl in Andromeda 
could be within the range of GLAST and ACT experi- 
ments if the DM cusp has steepened due to the influence 
of an IMBH. 

Gl suffers from a high uncertainty in the determination 
of its dark matter density profile; more measurements of 
velocity dispersion profile are needed, as well as a better 
understanding of the stellar mass-to-light ratio. Gl is 
also an extragalactic source, a factor of ~ 10 farther then 
typical dwarf galaxies of the Milky Way. On the positive 
sides, it has potentially an IMBH in its center which could 
substantially enhance the dark matter signal. Because 
of its old stellar population, Gl should not have a high 
astrophysical background in high energy gamma rays. 

With the GLAST detector coming up in the next year 
many new gamma ray sources will be discovered. It is 
possible that some of them will be uniquely interpreted 
as a dark matter signal. If we are lucky, the signal from 
Gl could add some information in this respect and maybe 
help disentangling the dark mattcr-astrophysical proper- 
tics. 
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